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Abstract. We study three-state Potts spins on a square lattice, in which all 
bonds are ferromagnetic along one of the lattice directions, and antiferromagnctic 
along the other. Numerical transfer-matrix are used, on infinite strips of width 
L sites, 4 < L < 14. Based on the analysis of the ratio of scaled mass gaps 
(inverse correlation lengths) and scaled domain-wall free energies, we provide 
strong evidence that a critical (Kosterlitz-Thouless) phase is present, whose 
upper limit is, in our best estimate, Tc = 0.29 it 0.01. From analysis of the 
(extremely anisotropic) nature of excitations below Tc, we argue that the critical 
phase extends all the way down to T = 0. While domain walls parallel to the 
ferromagnetic direction are soft for the whole extent of the critical phase, those 
along the antiferromagnctic direction seem to undergo a softening transition at a 
finite temperature. Assuming a bulk correlation length varying, for T > Tc, as 
5(T) = exp l^fej (T — Tc)~''j, cr ~ 1/2, we attempt finite-size scaling plots of 
our finite-width correlation lengths. Our best results are for Tc = 0.50 it 0.01. 
We propose a scenario in which such inconsistency is attributed to the extreme 
narrowness of the critical region. 



PACS numbers: 05.20.-y, 05.50.+q, 64.60.Fr, 75.10.Hk 

1. Introduction 

Many spin systems exhibit exotic low-temperature phases, frequently as a result of 
competing interactions coupled with suitable ground-state degeneracy. In this paper 
we study a borderline case in which, although neither frustration nor macroscopic 
ground-state entropy are present, a critical phase (i.e. with power-law decay of 
correlations against distance) arises for a finite temperature range above zero. The 
prototypical system in which such sort of phase occurs is the two-dimensional XY 
model The corresponding transition from the high-temperature (paramagnetic) 
to the critical phase, characterized by an exponential divergence of the correlation 
length, is known as Kosterlitz-Thouless (KT) transition. Although the continuous 
symmetry of XY spins plays an important role in the formation of vortices, whose 
presence is an essential feature of the low-temperature behaviour, it was soon realized 
that critical phases could be sustained also in some two-dimensional magnets with 
discrete spin symmetry. Examples are Z{N) models with N > 4 and 
models with a macroscopically degenerate ground state to which suitable degeneracy- 
lifting interactions are added. The latter case is illustrated by the triangular Ising 
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antiferromagnet in a uniform magnetic field ||, ^, or with ferromagnetic second- 
neighbour interactions @, ||, |^ • 

We consider the uniformly anisotropic Potts model with nearest-neighbour 
couplings on a square lattice pO| , whose Hamiltonian reads: 

H — - ^ ( Jx S{ax,y — ax+l,y) + JyS{ax,y - <^x,y + l) ) (1) 

x,y 

where the Potts spins (Jx^y take any one of q values, and the lattice site coordinates 
are given by {x,y). 

Our interest is in the mixed (ferro-antiferromagnetic) case ^ 12 J^; > 0, 
Jy < 0. For any q, the Hamiltonian is unfrustrated, as the product of interactions 
around plaquettes is always positive. The ground state of an L x L system has 
degeneracy 0{q- 1)^; thus, even though this diverges in the thermodynamic limit 
for q > 2, the residual entropy per spin approaches zero as L~^. 

The g = 3 mixed model which will concern us here has attracted special 
attention ^ as it is believed to exhibit a KT transition. 

Furthermore, evidence has been found that no transition occurs for q > 3 These 
two findings may in fact be viewed as complementary [T^ . In previous work, a great 
deal of effort was dedicated to calculating the critical temperature Tc above which the 
system is paramagnetic, particularly the variation of Tc against the ratio Jx/ Jy, with a 
rather large degree of scatter in the corresponding results |l4| [T5[ |l6[ [T7| . Although, 
in these latter references, assorted evidence is quoted in favour of the view that the 
transition is indeed "of an unconventional type" p^ , a more precise characterization 
is attempted only in , where the correlation lengths from Monte Carlo studies of 
L X L systems, L < 60, are shown to fit reasonably well to the exponential form in 
(T — Tc)~^^^, expected to hold close to a KT transition An investigation of the 
problem in its (l-l-l)-dimensional quantum version ||l3| lends numerical support to the 
respective form for the case, namely the mass gap vanishing as exp (—a (A — Ac)^°^), 
where A is the temperature-like coupling. However, the combined effects of small 
system sizes and uncertainty in the determination of Ac allow the authors of IT^ to 
conclude only that cr < 0.9. 

As regards the low-temperature region, the critical (or massless) phase is analysed 
in some detail in upon consideration of the Roomany-Wyld approximants [p^ 
to the beta-function jl9[ it is deemed as plausible that the mass gap vanishes 
all the way from A = to Ac (corresponding, in the classical analog, to a critical 
phase from T = to Tc). In a Monte Carlo rcnormalization group study of the 
chiral Potts model |2l] (with the parameters set in such a way as to map onto the 
present model the correlation functions are shown to behave consistently with 

the existence of a critical phase. However, the temperature range considered in [pT| 
goes only halfway down to zero from Tc- 

Thus, it seems worthwhile investigating the behaviour of this system in more 
detail, both (a) at the transition and (b) at low temperatures. In connection with (b), 
it must be recalled that, for triangular Ising antiferromagnets with second- neighbour 
ferromagnetic interactions 0, ||, ||, the critical phase does not extend all the way 
down to T = 0: there is an ordered phase with nonzero magnetization for a finite 
temperature range below the KT phase. While strictly similar behaviour should not 
be expected in the present case (since here the ground state has zero macroscopic 
magnetization), the question still remains of precisely how the (highly anisotropic) 
correlations spread at low temperatures. In pi], n2[, arguments are given to show 
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that slightly above T = an effective ferromagnetic coupling arises between second- 
neighbour rows of ferromagnetically coupled spins. This would single out a small 
subset of row configurations, namely 010101, 020202 and 121212 as energetically 
more favourable. It appears that the effect of this on the free energy (hence, on the 
corresponding thermodynamic state) has not been quantitatively investigated so far. 

We apply numerical transfer-matrix (TM) methods to the q ~ 3 mixed Potts 
model of equation (|l|), on strips of a square lattice of width 4 < L < 14 sites 
with periodic boundary conditions across (except where domain- wall free energies are 
analysed, in which case both periodic and twisted boundary conditions are considered). 
We shall restrict ourselves to Jx/Jy ~ — 1- 

In Section || we consider the scaling of correlation lengths, as obtained from the 
two largest eigenvalues of the TM. along both directions x and y of equation (|l|). From 
these we calculate the ratios of scaled gaps, which provide relevant information on the 
extent of the critical phase. In Section domain-wall free energies are calculated, and 
their associated correlation lengths are investigated; in Section ^ the region near the 
transition from critical to paramagnetic behaviour is examined, especially as regards 
the form in which the thermodynamic correlation length diverges as T — s- T+ ; finally, 
in Section |^, concluding remarks are made. 

2. Correlation lengths 

Here we consider the largest correlation length £,l{T) on a strip of width L, i.e. the 
one given by ^2^{T) = ln(Ao/| Ai |) where Aq.i are the two largest eigenvalues of the 
TM in absolute value. Thus it picks up the slowest decay of correlations along the 
specified direction of iteration of the TM. 

With the correlation function G(r) = ((5((To — Cr)) — 1/3 for the ferromagnetic 
direction one has G(r) ~ exp(— r/^^), while in the antiferromagnetic direction 
G(r) ~ (-l)-- exp(-r/e^^). 

As remarked in , the phenomenogical-renormalisation equation for the critical 
temperature, 

L^2\T*) = L'^Z'(Tn , (2) 

only admits a fixed point T* when the TM is taken along the ferromagnetic direction. 
The raw data for i^^^ against temperature are shown in Fig. respectively (a) for 
the TM taken in the ferro- (F) and (b) antiferromagnetic (AF) directions. In the 
latter case, curves for consecutive strip widths get almost tangent as L grows, in the 
region T ~ 0.5 — 0.6 (the L = 12 and 14 curves reach within one part in 10^ of each 
other, at T = 0.5) . A clearer understanding of both sets of data can be gained as 
follows. 

For the case of Fig. ^ (a) we plot the ratio of scaled mass gaps [|[ ^ ^ U| 

Qll'{T)^ [L^l\T)]/[L'^l}iT)] , (F) (3) 

where we always take L' = L + 2. At the fixed point of equation (||), Qll'{T*) = 1, 
and one expects it to remain close to this value for an extended temperature interval 
if a critical phase is present |l^. |2^, |2^, |2j]. Indeed this is what happens here, 
as shown in Fig. |2[ Starting slightly below the respective T*, and down to the 
lowest temperatures reached before running into numerical problems (typically T = 
0.03 — 0.04), the curves remain rather fiat. The values of Ql u at which they stabilize 
are, respectively, 1.04724, 1.01596, 1.00730, 1.00394, 1.00239, for L = 4, 6, ... 12. This 
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Figure 1. L^j^ against temperature: (a) TM in tlie ferromagnetic direction; 
(b) TM in the antiferromagnetic direction. Bottom to top; L = 4, 6, . . . 14. 



sequence behaves as Ql l' = a + bL~^, where x ~ 2.7 and a equals unity within one 
part in 10'^. As regards the extent of the flat region in the thermodynamic hmit, the 
shape of curves in Fig. ^ suggests that the upper hmit be given by extrapolation 
of the sequence of fixed points of equation (^): = VmiL^aoTI j^,. Such sequence is 
the same given in Table 1 of with the additional point Tj*2 14 — 0.40030. Using 
L = 4, 6, ... 12 we get T^^, ^T^ + c L^y with = 0.29 ± 0.01 (to be compared with 
Tc = 0.32 ± 0.03 0), y — 0.7. Thus we have convincing evidence that the system 
exhibits a critical phase, which extends down from Tc — 0.3, at least to the lowest 
temperatures examined, and presumably all the way down to T = 0. 

In order to check on this latter assumption, it is instructive to look at the low- 
temperature behaviour of the correlation lengths themselves. This is illustrated in 
curve (a) of Fig. ^ which shows that, immediately below Tc the functional form 

= f[l) exp(-l/T) (with F{L) - i^i, from the results for Qll'{T) above) sets 
in along the F direction. This corresponds to the lowest-energy (T ~ 0) excitations 
in that direction, which are obtained by breaking a single ferromagnetic bond, eg 
. . . 000000 . . . ^ . . . 000111 . . . (while, in this example, both neighbouring chains are 
. . . 222222 . . .). See below the contrasting case for the AF direction. In other words, 
T ^ behaviour is dominant all the way down from Tc, consistent with the idea that 
there is a single (critical) phase from T = upwards to T ~ 0.3. 

Examination of data taken with the TM along the AF direction. Fig. |^ (b), is 
helped because the zero-temperature bulk correlation length in that direction is known 
exactly 

r^(r = 0)=ln2. (AF) (4) 
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Figure 2. Qll'{T) of equation 
L = 4,6,. . . 12 



against T, for L' = L + 2. Top to bottom: 



We then look at the finite-temperature and finite-size deviations of (,2^iT) from 
equation (||), as displayed in curve (b) of Fig. ^. We have been able to reach roughly 
the same temperatures (usually T ^ 0.03) as with the TM along the F direction, 
before running into numerical problems. The trend shown for L = 14 in the Figure is 
followed for all strip widths, namely below T ~ 0.3 one has 

In 2 - Cl\T) = a{L) exp(-2/r) . (AF) (5) 

We have found empirically that the L— dependent prefactor is well fitted by a form 
a{L) = oq cxp(— ao ~ 500, h ~ 19.3, except for L = 4 which is off by ~ 50%. 

The factor of two in the exponential reflects the single-spin nature of elementary 
excitations along the AF direction, as the least costly move in this case consists in 
overturning a spin (breaking both its ferromagnetic bonds) without breaking any of 
its antiferromagnetic ones. Again, the fact that such low-T behaviour takes over 
immediately below T^. indicates that there is a single phase between T = and T^.. 



3. Domain- wall scaling 

For bulk d— dimensional systems, hyperscaling implies that close to the transition, the 
domain- wall free energy AF vanishes as the correlation length ^ diverges, according to 

id- - 



AF 

T - 



T: 



-1) 



|25 



Assuming power-law singularities, ^ - iTc-T]-", AF ~ (Tc-Ty, 
one has fi = (d — l)t^ ||2^. Specializing to d = 2 and using standard 
finite-size scaling ]2^ , one can derive the basic fixed-point equation of domain- wall 
renormalisation group on strips p7[ : 

LAfL{T*)=L'AfL'{T*), (6) 
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Figure 3. Curve (a): against 1/T, L = 14, TM in the F direction. Squares 
are calculated points; connecting curve approaches behaviour ~ exp(— 1/T) from 
1/T ~ 3 onwards. Curve (b): In 2 - ^"^ against 1/T, L = 14, TM in the AF 
direction. Triangles are calculated points; connecting curve approaches behaviour 
~ exp(-2/T) from 1/T ~ 3 onwards. 



where A/l is the free energy per unit length, in units of T, of a seam along the full 
length of the strip: A/l = fT^fh^ with (/J) being the corresponding free energy 
for a strip with periodic (twisted) boundary conditions across. 

While, for Ising spins, twisted and antifcrromagnetic boundary conditions 
coincide |2^, careful consideration ||2^ shows that in Potts systems with q > 2 this 
does not hold. Indeed, in order to produce a seam along the strip one needs an 
interaction which favours pairing between well-defined (and different) spin states, e.g. 
1 — 2, 2 — 3, 3 — 1. An antifcrromagnetic Potts interaction only makes selected pairings 
energetically unfavourable, while leaving all other combinations equally probable. 

Therefore, for Potts models with q > 2, twisted boundary conditions are obtained 
by changing the interaction across the L— direction as follows |28|: 

JS{ai-aL) — > J(5(o-i - (Tl + 1 (mod gr)) . (7) 

With the above definitions, one has A/l = — ln(AQ /A^) where Aq , are the largest 
eigenvalues of the TM, respectively with periodic and twisted boundary conditions 
across. In the present case, different results are to be expected, depending on whether 
the TM is set along the F or AF directions. In Figure ^ we show LA/^ against 
temperature for both cases. 

For case (a) , the low-temperature behaviour is similar to that found in Section H, 
in that all curves (i) very nearly collapse on top of each other, and (ii) behave as 
exp(— l/T). Analogously to the ratio of scaled mass gaps Qll' used in that Section, 
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Figure 4. L A/i, against temperature, (a): TM in the ferromagnetic direction. 
On this scale, curves are indistinguishable up to T ~ 0.4; top to bottom (on 
right-hand side): L = 4, 6, . . . 14. (b) TM in the antiferromagnetic direction. 
Curves cross in the central region of T axis. Top to bottom (on right-hand side): 
L = 4,6, ...14. 



we define the ratio of scaled free energies: 

nLL'{T)^[LAfL{T)]/[L'AfL'iT)], (F) (8) 

Results, using L' = L + 2, are displayed in Figure || and are very similar to those shown 
in Figure ||, except that instead of crossing the TZ — 1 axis from below, the curves only 
get asymptotically close to it, and from above. This is expected, as inverse correlation 
length and domain- wall free energy are dual to each other |^ . The values of TZl l' 
at which the curves stabihze are, respectively, 1.03057, 1.01029, 1.00470, 1.00253, 
1.00152, for L = 4, 6,... 12. This sequence behaves as TZll' = e + fL~'', where 
z ~ 2.7 and e equals unity within two parts in 10^. Though the absence of a crossing 
of the TZ = I axis deprives one of an obvious reference point, it is possible to infer the 
limiting extent of the flat region, as L oo, as follows. For each curve one records at 
what temperature it reaches within, say, 10~^ of its own T — > asymptotic value, and 
the resulting sequence is extrapolated against 1/L. By doing so we get Tc — 0.30 ±0.04 
where the error bar, though somewhat subjective, is certainly generous enough to allow 
for the degree of arbitrariness involved in this procedure. 

We conclude that the evidence provided is entirely consistent with that given by 
the analysis of correlation lengths in the F direction. As a final remark, recall that 
with the TM in the F direction, the seam whose energy is calculated lies parallel 
to ferromagnetic chains. With q = 3, the substitution in equation allows for a 
zero-energy seam at T = 0, hence the behaviour depicted in Fig. (a). 
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Figure 5. TZll'{T) of equation (jsj) against T, for L' = L + 2. Top to bottom: 
L = 4,6,. . . 12 



As regards data displayed in Figure Q (b), one finds well-defined fixed points 
of equation (||). With L' = L + 2, these are located respectively at T* = 0.58753, 
0.52589, 0.48334, 0.45108, 0.42543, for L = 4,6,...12. When plotted against 
this sequence displays accentuated downward curvature. Attempts to fit it against 
L~* produce t ~ 0.1, with a negative extrapolated value for T^. A two-point fit of 
L 10 and 12 data against gives an estimate = 0.29. Together with the 
downward curvature of the full set of data, this indicates that the actual given by 
this particular domain-wall scaling is below 0.29. 

We recall that in the previous case with the TM along the F direction, the 
corresponding domain walls were found to be critical for all < T < Tc ~ 0.30. With 
the TM along the AF direction the seam runs perpendicular to that least energetic 
orientation of domains, breaking ferromagnetic bonds as it goes. It would thus appear 
that, if there is a well-defined temperature for the softening of this latter type of 
domain walls, it should be finite. However, we have not managed to produce a reliable 
estimate for this from our sequence of fixed points. Bearing in mind that the only 
typical non-zero temperature arising so far in the problem is T,, ~ 0.30, one might 
interpret the data as signalling that the extrapolated critical point for domain walls 
along the AF direction is exactly T — 0. This brings new problems, as consideration 
of the physical domain- wall free energy per site, AF^ = TAf^/L shows that it 
behaves as 1 — a{L)T for T ^ 0, with a{L) ~ 2.7 exp(— 4.6/i). In other words, 
at low temperatures it becomes extremely expensive to excite a domain wall along 
this direction. While this behaviour is, in a loose sense, dual to that of the correlation 
length along the AF direction, which remains finite at T = as given by equation (^) , 
it is hard to reconcile with the idea that the respective domain walls would go soft at 
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T — 0. At present we are not able to offer a consistent interpretation for the beliaviour 
just exhibited. 

4. Approach to Tc from above: finite-size scaling 

Next we attempt a characterization of the transition from the paramagnetic to the 
critical phase. It is usually assumed that, as this point is approached from above in a 
Kosterlitz-Thouless transition, the bulk correlation length behaves as pO| : 

e(T) = exp (T - T,)-^] , (9) 

with a = 1/2, and a^, and Tc are non-universal. From finite-size scaling [26| , the 
correlation length on a strip is expected to scale as: 

LCZ'(r) = /(Lexp[-65(T-T,)--]) . (10) 

We have taken our data for correlation lengths calculated with the TM along the 
F direction, displayed in Fig. || (a) as LS^^^iT), and attempted to fit them to the 
form given by equation (|l0|). We used the procedure described in [Q, to try and 
collapse our results for assorted values of L and T onto a single curve. We took as 
basic sets of data the six sequences taken each for fixed L = 4, 6, ... 14 against varying 
temperature, from T = 0.04 to 1.1. According to ||3l|, we then took each of these 
sets in turn as a baseline, against which to fit the five other sets. Since equation ( p^ 
is valid only above Tc, the actual number of data used would vary, for each trial fit, 
depending on the particular value taken for the critical temperature. With Tc in the 
range 0.30 — 0.50, as explained below, the number of points fitted was between 280 
{Tc ~ 0.30) and 230 (Tc ~ 0.50). For the local fitting procedure fSljl, we used cubic 
polynomials. 

We started by keeping cr = 1/2 fixed, and varying and Tc- From Sections ^ 
and I Tc = 0.30 would seem a reasonable initial guess. Inspired in work done on the 
XY model fej was initially assumed of order unity. The landscape we found in 
parameter space, for the sum of residuals [ pl| as a function of &^ and Tc, was similar 
to that described in pO[ |, with "a number of local minima ... at which ... minimization 
programs can get stuck" . Furthermore, visual inspection of plots produced with the 
parameters set at the minimizing values gave consistently poor fits in all cases for 
which Tc was below 0.35. We did, however, manage to find an overall minimum for 
Tc = 0.51 ± 0.01, b^ = 1.72 ± 0.01 (with none for 0.35 < Tc < 0.50). 

We then allowed a to vary. The overall picture was entirely similar to that just 
described for fixed a = 1/2. The only solution which gave a visually confirmed good 
fit was Tc = 0.50 ± 0.01, b^ = 1.60 ± 0.01, cr ~ 0.54. The corresponding results are 
shown in Fig. pi The value Tc ~ 0.50 is inconsistent with the evidence exhibited in 
Sections ^ and]|. One possible solution to this contradiction would be if the critical 
region, where the behaviour of equation (|l^) is valid, were actually very narrow. Thus, 
we would be distorting the picture, by using high-temperature data which in fact are 
outside the critical region. For the XY model in two dimensions, it has been argued 
that this is indeed the case |3^, ^ , with the width of the critical region estimated 
as (T - Tc)/Tc < 10"^ Hi. With Tc ~ 0.3, in the present case this would mean a 
temperature interval AT ~ 3 x 10"^. For the fitting shown in Fig. |6| we used data 
for T as high as 2.5, for the narrowest strips. We noticed that, if only a subset of 
data with T < 1.1 is used, the best fit is found for Tc ~ 0.47, and is of similar quality 
to that of Fig. |^. Thus a plausible scenario is one in which, as only data for an ever 
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L exp [-b/(T-Tj''] 



Figure 6. Scaling plot (see equation (|lO|)) of data displayed in Fig.|l| (a). The 
parameters are: Tc = 0.51, 6j = 1.60, cr = 0.54 (see text). 



narrowing interval above are considered, the corresponding fits give rise to lower 
and lower estimates of itself. Should the prediction of Q hold true in the present 
case, one would have a problem in that the accuracy to which is known so far 
(taking, say, our best estimate 0.29 ± 0.01 of Section H) is lower than the width AT 
of the critical region. 

5. Conclusions 

We have provided strong numerical evidence that a critical (Kosterlitz-Thouless) 
phase is present in the three-state mixed Potts magnet described by equation (|l|). 
Our argument relies on the analysis of the quantities Qll'{T) and TZll'{T), defined 
respectively in equations (||) and (||). Based on data exhibited in Figures |^ and || and 
their respective extrapolation, we conclude that the upper limit of the critical phase, 
above which the system behaves paramagnetically, is Tc = 0.29 ± 0.01 (from Figure H; 
treatment of data in Figure ^ produces the less accurate, but not inconsistent, estimate 
Tc = 0.30 ±0.04). 

Analysis of the (extremely anisotropic) nature of excitations below Tc, as depicted 
in Figure ^, gives credence to the idea that the critical phase extends all the way down 
to T = 0, similarly to what happens in the two-dimensional XY model. 

We have found that, while domain walls parallel to the ferromagnetic direction 
are soft for the whole extent of the critical phase, those along the antiferromagnetic 
direction seem to undergo a softening transition at a well-defined, finite temperature. 
However, from our data we have not been able to produce a reliable estimate for 
location of such transition. 
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Attempts to fit our data for T > to the finite-size scaling form of equation (|o|) 
liave met with qualified success. While we have been able to produce rather good 
collapse plots, as depicted in Figure ^, this has been done at the price of setting 
Tc = 0.50 ± 0.01 which is clearly inconsistent with our own earlier estimates. We 
have proposed a scenario in which such result is interpreted as a distortion due to the 
narrowness of the critical region, similarly to the case of the XY model |3^, ^ Q . 
This is qualitatively borne out by the fact that fits with only a subset of data, taken 
at not too high temperatures, indeed result in a lower estimate for Tc. 
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